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We estimate the energy of the lowest charmonium and upsilon states with hybrid admixtures 
using the method of QCD Sum Rules. Our results show that the ^'(2S) and T(3S) states both 
have about a 50% admixture of hybrid and meson components. From this we find explanations of 
both the famous p — n puzzle for charmonium, and the unusual pattern of a decays that have been 
found in T decays. Moreover, this picture can be used for predictions of heavy quark production 
with the octet model for RHIC. 



PACS Indices:14.40.Gx,12.38.Aw,11.55.Hx,13.25.Gv hybrid states for the RHIC search for the QCD phase 

transition via heavy quark state production. In Sec. 
V we review our conclusions. 



I. INTRODUCTION 

There is a great interest in studying states with 
active glue, such as hybrid mesons, a color singlet 
composed of a quark-antiquark in a color octet and 
a gluon, in order to better understand nonpertur- 
bative QCD. Recently we have used the method of 
QCD Sum Rules in an attempt to find the lowest hy- 
brid charmonium state pj. Our conclusion was that 
the physical states with active glue must be mixed 
states, with both charmonium and hybrid charmo- 
nium components. In the present work we use QCD 
Sum Rules for J PC = 1 vector states to find the 
lowest mixed meson-hybrid meson states for both 
charmonium and upsilon systems. 

In addition to the importance of finding states 
with active glue, we are motivated by several exper- 
imental considerations. First, the ratio of hadronic 
decays of the charmonium ^f'(2S) compared to the 
J/ , 5(1S') state is more than an order of magni- 
tude smaller than predicted by perturbative QCD 
(PQCD), the so-called p — tt puzzle, which was 
discussed at length in Ref[l[. Second, the T(nS) 
states have an unusual pattern of decays into two pi- 
ons, which also cannot be consistent with PQCD [2], 
which we call the Vogel T(An = 2) puzzle. Thirdly, 
our theory of heavy quark states provides a basis 
for the color octet model predictions of RHIC heavy 
quark production^, 0| 

In Sec. II we discuss the motivation for the present 
work on mixed heavy quark and hybrid mesons. In 
Sec III we review the method of QCD sum rules, the 
work in Ref[l[ on hybrid charmonium, and apply the 
method of QCD Sum Rules for mixed meson-hybrid 
meson charmonium and upsilon states. In Sec. IV 
we discuss our solution to the p — it puzzle, the Vogel 
T(An = 2) puzzle, and applications of our mixed 



II. HEAVY QUARK PUZZLES AND RHIC 
EXPERIMENTS 

First let us look at the lowest charmonium and 
upsilon (nS) states (FIG 1): 

Y(4S) 10,580 

Y(3S) 10,355 

▼PS) - - 3686 Y(2S) 10 023 



J/V(1S) 3097 Y(1S) 9460 

cc states b ^ states 

FIG. 1: Lowest nS charmonium and upsilon states 



Note that the separation in energy between the 
ip'(2S) and J/ip(lS) states is nearly the same as 
the separation energy of the T(2S) and the T(IS') 
states. This will be important for our studies of 
Heavy Quark Hybrids, but turns out to be mislead- 
ing. 



A. The p — n Puzzle 

The p — 7T puzzle for cc 1 states is based on the 
two diagrams for PQCD and electromagnetic decay 
of such states, shown in Fig. 2. By taking ratios the 
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FIG. 2: Perturbative QCD and em diagrams for 
QQ(1~~) Decays 



wave functions at the origin cancel, and this pre- 
dicts the ratio of branching ratios for cc decays into 
hadrons (h) 



B(V'(2S) -> h) B(^'(2S) -* e+er) 

£(J/tf (IS) -> h) ~ B(J/tf(lS) -> e+e-) 
~ 0.12, (1) 

the famous 12% rule. 

The p - 7T puzzle: The W(2S) to J/* ratios for 
p — 7r and other hadron decays are more than an 
order of magnitude smaller than predicted by Eq|lj). 
Many theorists have tried and failed to explain this 
puzzle. See Ref|5] for a review. This and more recent 
attempts at a solution are also discussed in RefQ, 
and all agree that previous work has not produced 
a solution for this puzzle. 



C. The Octet Model for RHIC and Hybrids 

The major goal of modern RHIC (Relativistic 
Heavy Ion Collision) experiments is to produce and 
study the quark-gluon plasma (QGP) which existed 
in the early universe before the QCD phase trasition, 
about 10 -5 seconds after the big bang. One impor- 
tant signal of this QGP is the production of heavy 
quark (charmonium and upsilon) states via qq in- 
teractions in the early universe. The most natural 
mechanism is qq — > g — > QQ, in which an octet 
qq produces an octet QQ, which is just PQCD, fol- 
lowed by the nonperturbative (NPQCD) process in 
which the octet QQ becomes a singlet QQ with the 
emission of a gluon (or other color octet). This is 
depicted in Fig. 3: 




< ---- Perturbative QCD »--> <--- Nonperturbative P ----> 



FIG. 3: Octet model for production of heavy quark 
mesons 



B. The sigma Decays of T(nS) States Puzzle 

The puzzle of sigma decays of bb 1~~ (T(nS)) 
states is given by the following. The sigma is a 
low-energy broad two-pion scalar resonance. Exper- 
iments on Tfn S) states find0 

T(25) — > T(15) + 2ir has a large branching ratio, 
but no a 

T(3s) — > T(IS') 4- 2ir has a large branching ratio 
to a 

An = 2, emit a 

An 2, no a emitted. 

This is the Vogel An = 2 puzzle, which cannot 
be understood using perturbative QCD, as expected 
for heavy bottomium states. 



The nonperturbative matrix elements for the tran- 
sition from the color octet < QQ(8) \ state to a color 
singlet ^ state, < 1 1 > in the notation of Ref[6], 
have been determined by fits to experiments using 
the octet modelQ. As we shall see, our determina- 
tion of mixed heavy quark and heavy quark hybrid 
mesons will provide a mechanism for predicting these 
NPQCD matrix elements. 



III. MIXED HEAVY QUARK HYBRID 
HEAVY QUARK 1~ STATES AND QCD 
SUM RULES 

In this section we review the method of QCD sum 
rules, review our previous application of this method 
to attempt to find the lowest energy hybrid charmo- 
nium 1 state, and present our new application of 
the QCD sum rule method to find the lowest energy 
mixed charmonium and upsilon states with hybrids. 



3 



A. Method of QCD Sum Rules 

The starting point of the method of QCD sum 
rules 0] for finding the mass of a state A is the cor- 
relator, 



Il A (x) = (\T[J A (x)J A (0)]\) 



(2) 



with |) the vacuum state and the current Ja{x) cre- 
ating the states with quantum numbers A: 



Ja{x)\) = c A \A) +J2 c n\n; A) , 



(3) 



where \A) is the lowest energy state with quantum 
numbers A, and the states are higher energy 

states with the A quantum numbers, which we refer 
to as the continuum. 

The QCD sum rule is obtained by evaluating H A 
in two ways. First, after a Fourier transform to mo- 
mentum space, a dispersion relation gives the left- 
hand side (lhs) of the sum rule: 



n(<?)i£ s 



Imn A (M A ) 



q 2 ) 



7r(s — q z ) 



where Ma is the mass of the state A (assuming zero 
width) and s a is the start of the continuum a pa- 
rameter to be determined. The imaginary part of 
H A (s), with the term for the state we are seeking 
shown as a pole (corresponding to a S(s — M\) term 
in ImlT) , and the higher- lying states produced by J a 
shown as the continuum, is illustrated in Fig. 4: 



ImflKs 




fil to continuum 



FIG. 4: QCD sum rule study of a state A with mass Ma 
(no width) 



Next H A (q) is evaluated by an operator product 
expansion (O.P.E.), giving the right-hand side (rhs) 
of the sum rule 



n(g) r l = 5> fc (g)<o|0 fc |o) , 



(5) 



where Ck(q) are the Wilson coefficients and (0|Cjt |0) 
are gauge invariant operators constructed from 
quark and gluon fields, with increasing k correspond- 
ing to increasing dimension of Ok ■ It is important to 
note that the Wilson coefficients, Ck(q) obey renor- 
malization group equations [8| 

After a Borel transform, B, in which the q variable 
is replaced by the Borel mass, Mb, 



B 



2 lim T^w^t-A: 

ij 2 ,n— >oo (n — 1)! dq^ 



(6) 



q*/n=M* 



the final QCD sum rule, BU A (q)(LHS) 
BU A (q)(RHS), has the form 



l e -Ml/M% +B f IjL 



Im[U A (s) 



■ds 



r(s - q 2 ) 

= BY / c A (q)<0\O k \0> . (7) 

k 

This sum rule and tricks are used to find Ma, 
which should vary little with Mb- A gap between 
M\ and s is needed for accuracy. If the gap is too 
large, the solution is unphysical, which is important 
for our present work, as we discuss below. 

B. Hybrid charmonium 

Here we give a brief review of the calculation of the 
correlator and the results for the QCD sum rule for 
a pure hybrid charmonium 1 state, which could 
possibly be the ^'(25). The current Jhh (which we 
called Jh in Ref[l|) for a heavy quark hybrid meson 
with J PC = 1— is 



(8) 



where ^> is the heavy quark field, T u — C^j v , 7„ is 
the usual Dirac matrix, C is the charge conjugation 
operator, and the gluon color field is 



£-~> 9 a 



(9) 



with A the SU(3) generator (Tr[X a X b ] = 25 ab ). The 
correlator 

ngnix) = <o\t[J^)Jhh(o)]\o>, (io) 

after a Fourier transform, was evaluated using the 
leading two operators in the operator product expan- 
sion, shown in Figs. 5 and 6. The scalar correlator 
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FIG. 5: Lowest-order term in Sum Rule 




FIG. 6: Gluon condensate term in Sum Rule 



n s is defined by W v {p) = [p^Pu/p 2 ~ g^)Ii v {p) + 
(p^/p 2 )U s (p). 

The leading term in the OPE for U. HII , Uf HH , 
corresponds to the diagram in Fig. 5. It is quite 
complicated, and the result is given in Ref0. Af- 
ter the Borel transform, II (Mb) is given in the 
Appendix, Eq(??). 

The second order term, corresponding to the op- 
erator with a gluon condensate shown in Fig. 6, 
has a scalar part JI^ hh (Mb) given in the Appendix, 
Eq(??). 

In Ref[l[ a solution to the QCD Sum rule was 
found for a charmonium hybrid state at the mass 
of 4 ,/ (25 l ), from which one would at first conclude 
that the ^'(25) is a pure hybrid 1 meson. How- 
ever, in order to satisfy the critrion that the solution 
is almost independent of the Borel mass a value of 
s = 60.0Geu 2 was needed. This would imply that 
the next excited state was 7 to 8 GeV, which is not 
consistent with the first state at only 3.66 GeV. Note 
that lattice QCD calculations found the first char- 
monium hybrid at about one GeV higher than our 
solution^, [nj, which is also consistent the ^f'(2S) 
not being a pure hybrid. 

This result, as well as the heavy quark puzzles 
and RHIC experiments discussed in Sec II, were the 
main motivation for the present work, in which we 
seek a solution for a mixed charmonium and hybrid 
charmonium state. 



C. Mixed charmonium-Hybrid charmonium 
States 

Recognizing that there is strong mixing between a 
heavy quark meson and a hybrid heavy quark meson 
with the same quantum numbers (as shown below), 
and that the fact that our pure hybrid charmonium 
solution was not a physical state, we now attempt to 
find the lowest J PC = 1 charmonium state with 
a sizable admixture of a charmonium meson and a 
hybrid charmonium meson. An appropriate mixed 
vector (J PC = 1 ) charmonium, hybrid charmo- 
nium current to use in QCD Sum Rules is 



J" 



bJ 



with 



H 



TV 

If — 



HH 



(ii) 



(12) 



where J H is the standard current for a 1 char- 
monium state, and J HH 1S t Qe heavy charmonium 
hybrid current given above in Eqs([8l[9]). 

Therefore the correlator for the mixed state: 



n%- HH (x) = <o\T[j^(x)r 



0> (13) 



IS 



^H-HH ( x ) 



n£W(s) 



(,;) 



b^ HH (x) 



HH V 
HHH 



(,:) 



+2&V1 
<0\T{j£(x)J H (0)}\0> 

<o\T[j HH (x)r HH (o)]\o> 

<Q\T[j£(x)r HH (0)}\0> . 



(14) 



For heavy quarks the gluon condensate is propor- 
tional to the quark condensate, and the renormal- 
ization group equations for the Wilson coefficients of 
the operator product expansions of Hh, Hhh, and 
U-hhh are similar for the terms considered hereQ. 

The heavy hybrid correlator, n^ J -(g 2 ) was pre- 
sented in the previous section. The operator product 
expansion for the standard heavy quark correlator, 
Ilff ' (<1 2 ) (RHS) is given by the diagrams shown in 
Fig. 7. 

The leading term for the quark correlator in mo- 
mentum space, with Mc the charm quark mass, is 



S(k) 



1 5 



2 
fJv 



d 4 k 

(2tt)4 
M c 



Tr[S(knS(p-k)Tf] 



fc 2 

7^75 



Ml 



(15) 
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FIG. 8: Meson- hybrid meson lowest order diagram 




Tr[G" A (0)G ,t4 (0)] . 



+ higher order terms 

Quark propagator 
— ^ 0— Quark condensate <qq > 

Gluon condensate < G G > 
FIG. 7: Heavy quark meson diagrams 



Noting that the charmonium quark condensate is 
very small, and that the gluon condensate term and 
all higher-dimensional terms are also small, 11^ (p) 
dominates the heavy quark correlator, TLg (p). Car- 
rying out the momentum integral in Eq (|15[) , and ex- 
tracting the scalar correlator we find 



n|(p) 



da- 



Qp 4 - 2ip 2 Ml 



(47r) 2 J (a - a 2 )p 2 - 
Carrying out the Borel transform we find 

+3K 2 {2z)] , 



Ml 



(16) 



n|(M B ) = ^Mlexp- 2z [^-K {2z) 



2^i(2^) 
(17) 



with z = Ml/Ml. 

Finally, for the ^-hhh term, the dominant dia- 
gram is shown in Fig. 8, in which the gluon from 
the Jhh operator is coupled to a quark, leading to 
the Jh operator. This is essentially the perturbative 
plus nonperturbative H-HH matrix element without 
condensates. 

Using the external field method, the leading term 
of Hj^ HH , corresponding to Fig. 8, is 



u hhhi(p) 



■at 



-Tr\ 



[a KS ,(f(+Mc)]- 



4 J (2tt) 4 l (k 2 - 
Cjxip 1 - ^ + M c )(C ltl ) T 



M 2 ) 



(p - kf 



M 2 C 



(18) 



After a Borel transform and extracting the scalar 
component of n^ HH1 , one finds that 



n 



HHH 



7r 2 n^(A/ B ) 



(19) 



Therefore thr right-hand side of our scalar corre- 
lator is 



IT 



H-HH 



(M B ) rhs 



(b 2 + 2TT 2 bVl.-b 2 )U s H (M B ) 
+ (l.-b 2 )U s HH {M B ). (20) 



The left hand side of the sum rule has the usual form 



(see Eq(7)) 



n 



H-HH 



(M B )ihs = Fe- M "-™! M * 



(21) 



(K0- 



K1M 2 B 



K2M% + KZM%~ 



with s , K0, Kl, K2, K3 parameters used to fit the 
continuum. Note that the meson and hybrid me- 
son states associated with the H and HH opera- 
tors are normalized independently, and the operators 
have different dimensions. We renormalize by calcu- 
lating NHH = J dM B U s H (M B )/ J dM B U HH (M B ). 
Henceforth, for U HH (M B ) we use NHH x 

n| ff (M B ) 

As in Ref |1| , we obtain the expression for the mass 
of the mixed heavy meson-hybrid heavy meson by 
taking the ratio of the derivative of the sum rule 
with res[ect to 1/M B to the sum rule, giving 

M H -hh = {[s o (K0 + K1M% + K2M% + K3M B ) 

+KIM$ + 2K2M% + 3K3M B ]e~^B 
+d 1/AP n s H _ HH } x{(K0 + KIM 2 + K2M% 



+K3M%)e^ - U s Hh h y 1 



(22) 



A key parameter in our numerical fits is the value 
of b. The solution for b=-0.7 was most successful in 
fitting the criteria for finding the mixed hybrid state 
using QCD sum rules. The range of b for which 
a satisfactory solution is obtained is b=— 0.7 ± 0.1, 
with the result for b=-.7 shown in Fig. 9. 



G 



13.64 I 1 1 1 

13.50 13.55 13.60 13.65 



13.70 13.75 



FIG. 9: Mixed charmonium-hybrid charmonium mass 
3.69GeV 



continuum, and the mixing parameter b is also cho- 
sen to give a solution in which mixed upsilon state 
mass is almost independent of the Borel mass. The 
result is shown in Fig. 10. 
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We find the mass of the lowest-energy mixed 
charmonium-hybrid charmonium to be about the 
energy of the V(2S) state, 3.69 GeV, with s o =20 
GeV 2 , b= -0.7 => 50-50 per cent charmonium-hybrid 
charmonium. It satisfies the criteria for about a 
fifteen per cent accuracy. The values of the other 
parameters are KO = -15.9, Kl = 0.224, K2 = 
-0.00015, K3 = 0.00009. The only solutions sat- 
isfying the sum rule criteria are those with the value 
of b about — .7 ± .1, so that we find the state to be 
about a 50-50 per cent meson-hybrid meson. As we 
shall see, this gives a solution to the p — tt puzzle. 



D. Mixed upsilon- hybrid upsilon states 

The calculation of the mixed upsilon-Hybrid up- 
silon meson mass is the same as that of the mixed 
charmonium-hybrid charmonium mass using QCD 
sum rules, except the charm quark mass (which we 
took as Mq = 1.8 Gev 2 ) is replaced by the bottom 
quark mass (which we take as M 2 = 25.0 Gev 2 ). 
In fact, the QCD sum rule method is more accu- 
rate for the calculation of upsilon states, since the 
bottom quark condensate is much smaller that the 
charm quark condensate, and the operator product 
expansion converges faster. 

Since we found that the *B'(2S) is a mixed char- 
monium state (see previous subsection) and as we 
noted earlier the separation in energy between the 
\p'(2S) and J/ip(lS) states is nearly the same as the 
separation energy of the T(2S) and the T(1S) states 
(see Fig 1), we would expect that the T(25) is a 50- 
50 mixture of upsilon and hybrid upsilon. This in 
not our result, as we shall now see. 

From the QCD sum rule one obtains the expres- 
sion given in Eq. (l22j) . except the charm quark mass 
is replaced by the bottom quark mass in the expres- 
sions for the right-hand side of the correlator. The 
parameters s a , K0, Kl, K2, K3 are chosen to fit the 



FIG. 10: Mixed upsilon-hybrid upsilon mass = 10.4 GeV 



We find the energy of the lowest mixed upsilon 
meson and hybrid upsilon meson state to be at 10.4 
GeV, approximately the energy of the T(3S) state 
(see Fig. 1.). The parameters are s — 120 GeV 2 , 
K0 = -50, 000., Kl = 70500., K2 = -605., K3 = 
—0.5165, and b~ -0.7 for a good solution. Thus we 
predict that the T(3S) state is a 50-50 percent ad- 
mixture of an upsilon meson and a hybrid upsilon 
meson. As we shall now see, from this we have ob- 
tained a solution to the Vogel An = 2 puzzle. 



IV. MIXED MESON-HYBRID MESON, 
HEAVY QYARK DECAY PUZZLES, AND 
OCTET MODEL 

In this section we show that the solutions for 
the mixed nature of the charmonium ip'(2S) and 
bottomonium T(35) states provide explanations for 
the decay puzzles and a basis for the calculation of 
the nonperturbative matrix elements needed for the 
octet model used in RHIC calculations. 



A. The p — tv Puzzle 

First note that the matrix element 

< Trp\0\ip'(ce, 2S) > for p - n decay of \cc(2S) > is 
given by the PQCD diagram shown in Fig. 11 

Next, the hybrid decay matrix element 

< Trp\0'\ip'(ceg,2S) > is given by the PQCD dia- 
gram shown in Fig. 12. 

As one can see from the diagrams, these matrix el- 
ements are almost equal in magnitude. Since we find 
that \¥'(2S) >~ -0.7\cc(2S) > +0.7\ccg{2S) >, so 
the charmonium and hybrid charmonium approxi- 
mately cancel, we obtain for all 2 hadron decays, 




including p + tt decay, 

B(V'(2S) ^ p + tt) 



R 



« 0.12 , (23) 



B{J/*(1S) -> p + rr) 
which is our proposed solution to the p — ir puzzle. 

B. a Decays of T(nS) States Puzzle 

The solution to the Vogel An = 2 puzzle is 
based on the aplication of the glueball/sigma model, 
based on the study of scalar mesons and scalar 
glueballspj], G3, which was motivated by the BES 
analysis of glueball decay [Hj], and our solution for 
the lowest mixed state to be the T(3S I ) state. 
The glueball/sigma model has been used for pre- 
diction of sigma production from glue created in 
hadron-hadron collisions[14] and the decay of hybrid 
barvons[la|. which is closely related to the puzzle of 
sigma decays from upsilon states. The key is the 
glueball-meson coupling theorem(l6j 



dxT[J G (x)J m (0)} 



32 
~9 



<qq>, (24) 



where < qq >= Quark Condensate, which is de- 
picted in Fig. 13. 

From this one can calculate the matrix element for 
sigma decay from a hybrid meson, using the diagram 
shown in Fig. 14. 



Just as scalar glueballs, such as the f o (1500), de- 
cay mainly into sigmas, the hybrid component of the 
T(3S) has a strong a decay branch, while we predict 
that the T(2S) two-pion decay to the T(1S) would 
have a very small a decay branch. Therefore, our 
solution for the T(3S) to be a mixed bb-bbg provides 
a solution to the Vogel An = 2 puzzle. Since our 
states are not normalized we cannot calculate the 
numerical value of the cross section, a subject for 
future research. 



C. Octet Model for RHIC 

As discussed in section II, the octet model, de- 
picted in Figure 3., is the dominant mechanism for 
production of heavy quark states from a quark-gluon 
plasma produced via RHIC. Let us consider the colli- 
sion of a nucleus A, e.g., lead or gold nucleus, with a 
similar nucleus. The differential cross section for the 
production of a charmonium state in a A-A collision 
in the color octet model is 

-^-[pp -► V(cc)] = J f q/A f q/A ^[qq - CC(8) 



ip(cc)} 



da 
~dt 



[qq -> CC(8) -> ip(cc)] = [perturbative QCD] 



x < 0\ot Q \0 > 



(25) 
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where f q /A is the momentum fraction carried by a 

quark in the nucleus A, and < 0|Og Q |0 > is the 
NPQCD color octet matrix element. In previous ap- 
plications of the model, the nonperturbative octet- 
singlet matrix element was taken from fits to other 
experiments 0, Q|. We, however, can determine the 
NPQCD matrix elements for an octet quarkonium 
pair to emit a gluon (octet) and leave a physical sin- 
glet quarkonium state, which is given by lis,].. This 
is shown in Fig. 15. 




color octet-singlet matrix element can be obtained. 
This nonperturbative matrix element can be used 
for studies of the production of heavy quark states 
in RHIC experiments, using the octet model. It 
can also be used with the sigma/glueball model to 
predict the cross sections for sigma production fron 
heavy quark state decays. Since the states used in 
the QCD sum rule method are not normalized, these 
numerical estimates cannot be made at the present 
time 

In the near future we plan to extend our calcula- 
tion, so that numerical predictions of heavy quark 
decays and RHIC production of heavy quark states 
can be made. This will include possible tests of 
RHIC quarkonium production via sigma decays. 



VI. APPENDIX 



FIG. 15: Hhhh = ng,i oc color Octet-color singlet ma- 
trix element 



In RefQ we found that Uf HH (M B ), the Borel 
transform of the scalar term of the the main dia- 
gram for the HH correlator, shown in Fig. 5, is 



We have seen how to evaluate this diagram, but 
the normalization of the states must be carried out 
to make a numerical estimate. We can, however, 
estimate ratios of matrix elements, to predict ratios 
of quarkonium production. As an example, from a 
table in Cho-Leibovich[6J 

< 0|Og /V, (15)|0 >= 1.2 x lQ- 2 GeV 3 



< 0\Of (25)|0 >= 0.73 x lQ- 2 GeV 3 or 



Rc 



< 0\of(2S)\0 > 
<O|08 /V, (1S)|O> 



0.6 



(26) 



Since in our model of the ifj'(2S) state the cc(l) 
component should dominate, the parameter b~ — .7 
gives a rough estimate of this ratio, in agreement 
with the Cho-Leibovich phenomenological fit. 



V. CONCLUSIONS 

In summary, we find that the ip'(2S) is approxi- 
mately 50% charmonium and 50% hybrid charmo- 
nium; and the T(3S) is approximately 50% bot- 
tomium and 50% hybrid bottomium. This solves 
the p — 7r problem for charmonium decays, and the 
Vogel An = 2 puzzle for sigma decays of upsilon 
states. 

From the correlator corresponding to the mixed 
heavy meson and heavy hybrid meson current, the 



nWM B ) 



-gl 



2(4+) 



rM, 



{[-310 



. °° - 2 — 9- (1+5) 
4 ' dSe ^ 



638(5- 656(1 + 5)] 



1 + 5 

*s(2^|(l + 5)) + [-1892^ 



6065 - 1968(1 + <5) - 32(4(5 - 3 



+ 
S 2 



1 



6 3 



3(1 + 6) 



-4778- 



r 2 B 



-128(4(5-3 



M 2 



+ 94425- 4920(1 + 5) 
5 2 <5 3 „ 



1 + 5 3(l + 5) 2 



*i (2-^(1 + 5)) + [-1356 
+62845- 3280(1 + 5) 
-96(45 - 3— + 



1 + 5 



M 2 



^0(2^(1 + 5))} 
+multiple integrals 



(27) 



The multiple integral terms in Eq.(??) are small 
and are dropped. Mq is the charm quark mass for 



the charmonium calculations and the bottom quark 
mass for the upsilon calculations. We take Mq = 1.8 
GeV 2 and M h 2 = 25.0 GeV 2 . The gluon condensate 
term is shown in Fig. 6. After the Borel transform 
the scalar part of this term[l[, X\^ hh (Jv1b\ is 

M Q 

^Ihh(Mb) = - l 9 2 v^2 M Q e ^ ( 28 ) 
Mq 14 Ml Ml 

The K n are Bessel functions of imaginary argu- 
ment, related to Hankel functions by 

K n (x) = ^e-^H^iix) (29) 
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